Draft version February 2, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. ^jllj^A, 



GIANT MOLECULAR CLOUDS IN M64 

E. ROSOLOWSKY AND L. BlITZ 
Radio Astronomy Laboratory, University of California, Berkeley, CA 94720 
Draft version February 2, 2008 

ABSTRACT 

We investigate the properties of Giant Molecular Clouds (GMCs) in the molecule-rich galaxy M64 
(NGC 4826). In M64, the mean surface density of molecular gas is 2A^(H2) ~ 10^^ over a 2 kpc region, 
equal to the surface densities of individual GMCs in the Milky Way. We observed the J — 1 ^ Q 
transitions of CO, "CO, and HON. The line ratio Wco/W^i3CO for 200 pc < Rgai < 800 pc is 
comparable to that found in the Milky Way and increases significantly outside this region, in part due 
to a large contribution to the CO emission from diffuse gas, which composes 25% of the molecular mass 
in the galaxy. We developed a modified CLUMPFIND algorithm to decompose the ^'^CO emission into 
25 resolved clouds. The clouds have a luminosity-linewidth relationship L oc AU^'^^°'^, substantially 
different from the Milky Way trend reported by Solomon et al.l (^87): L oc AV^. Similarly, the 
clouds have a linewidth-size relationship of AV cx compared to AV cx in the Milky 

Way. Estimates of the kinetic and binding energies of the clouds suggest that the clouds are self- 
gravitating and significantly overpressured with respect to the remainder of the ISM in M64. The 
^■^CO-to-H2 conversion factor is comparable to what is seen in the Galaxy. The M64 clouds have a 
mean surface density at least 2.5 times larger than observed in Local Group GMCs, and the surface 
density is not independent of mass as it is in the Local Group: Sh2 cx m'^-'^^^-'^ . The clouds are 
correlated with the recombination emission from the galaxy, implying that they are star forming; the 
rate is comparable to that in other galaxies despite the increased densities of the clouds. The gas- 
to-dust ratio is similar to the Galactic value but the low extinction in the visual band requires that 
the molecular gas be clumpy on small scales. We note that the internal pressures of clouds in several 
galaxies scales with the external pressure exerted on the clouds by the ambient ISM: Pint oc Pg^p^*' ^^. 
We attribute the differences between M64 molecular cloud and those in the Local Group to the high 
ambient pressures and large molecular gas content found in M64. 

Subject headings: Galaxies:individual (M64) — galaxies:ISM — ISM:clouds — ISM:structure — radio 
lines:ISM 



1. INTRODUCTION 

Over 80% molecular gas in the Milky Way is found in 
Giant Molecular Clouds (GMCs) with masses 10^ Mq < 
Mgmc < lO^ *^ Mq and sizes between 2 and 100 pc. 
They have a colum n density of 2N(ll-?,) ^ 10^^ cm^^, in- 
dependent of mass llBlitdllQQ^ . Observations of molec- 
ular gas across the Local Group have found that most of 
the molecular gas in these systems is also found in GMCs 
with properties similar to those in the disk of the Milky 
Way iMizuno et al.ll200H lEngargiola et al.ll200'3t lYound 
l200nHBolatto et al.ll2nol . Beyond the Local Group, ob- 
servers find several galaxies where the typical molecular 
column density thro ugh the galactic disks is > 10^^ cm~^ 
ijHelfer et al.ll2003|l on kiloparsec scales implying that 
the ISM is overwhelmingly molecular, since the column 
densities of atomic gas are only rarely seen to exceed 
10^^ cm~^. Since this is larger than the typical column 
density through a single CMC in the Local Group, the 
molecular ISM must be qualitatively different in these 
galaxies. Do the molecular clouds in these galaxies blend 
together into a continuous sea of molecular gas or are 
there discrete clouds of molecular gas that are denser 
than Local Group GMCs? Moreover, what are the star 
forming structures in this molecular gas? In the Local 
Group, star formation is predominantly found in GMCs, 
but if there are no discrete clouds in this molecular gas. 
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do stars form uniformly throughout the galaxy from stel- 
lar mass clumps in the sea of gas? If there are GMCs but 
they are smaller and denser than Local Group clouds, is 
the star formation efhciency enhanced? In order to an- 
swer these questions, we have conducted a study of the 
molecule-rich galaxy M64 (NGC 4826) at high resolution 
to search for CMC analogues in a molecule-rich environ- 
ment. 

M64 is a good target for such a study because it has 
a mean column density across the central 2 kpc of the 
galaxy of 2A^(H2) = 10^^ cm~^, making it an intermedi- 
ate case between the Local Group and starburst galaxies. 
The galaxy is one of the closest molecule rich galaxies 
(4.1 Mpc, TuUv 1988). A merger with a small, gas-rich 
galaxy roughly 1 Gyr ago is thought to be responsible 
for the relatively l arge molecular g as content for its early 
Hubble type (SA) l|Braun et al.ll99 4. B94). The 3'.'5 res- 
olution obtained with the BIMA millimeter interferome- 
ter projects to 75 pc, which would be sufficient to resolve 
the largest molecular clouds in the Milky Way disk at the 
distance of M64. Two other recent studies examin e prop- 
erties of molecular emission in the nucleus of M64. iMeieil 
()2002l M02) conducted a study of the nuclear region us- 
ing the OVRO millimeter array in the 3 mm transitions 
of CO, i^CO, and C^^O. The primary beam at OVRO 
has a FWHM of 65" which samples gas to Rgai = 600 pc 
along the major axis and the map has a 5" synthesized 
beam. These observations were principally used to deter- 



2 



Rosolowsky & Blitz 



mine the bulk properties of the molecular gas using LVG 
analyses. iGarcia-Burillo et al.l (|2003, hereafter NUGA) 
also conducted an extensive analysis of the molecular 
gas in the nuclear region using the IRAM Plateau de 
Bure interferometer (PdBI) in the C0(2 1) and the 
C0(1 0) transitions. With the 42" primary beam, 
the PdBI observations cover the galaxy to Rgai = 420 pc 
for the (1^0) transition with a 2.5" resolution. The 
NUGA project focused interpretation of their data on the 
mechanisms by which the central AGN could be fed. Our 
search for molecular clouds in M64 is complementary to 
the aims of both M02 and NUGA. With a 105" primary 
beam and 3.5", our BIMA study maps the galaxy to the 
edge of the molecular disk at Rgai ^ 1000 pc permitting 
a more thorough study of the molecular gas. In addi- 
tion, this study includes zero-spacing observations of CO 
and ^•^CO which are essential in searching for low surface 
brightness gas among the clouds. In this paper, we re- 
port on our observations (fJSJi examine the dynamics of 
the molecular disk f H3.1(l and variations in the isotopic 
line ratios f ^3.2|l . decompose the blended emission into 
cloud candidates and argue that these candidates repre- 
sent the high density analogues of Local Group GMCs 

2. OBSERVATIONS 
2.1. Molecular Line Observations 

We observed M64 using the BIMA array l|Welch et al.l 
Il996|) in three tracers of molecular gas: C0(1 0), 
i3C0(l ^ 0) and HCN(1 ^ 0). We utilized data from 
the B, C and D arrays, including medium-resolution 
(C array) and zero-spacing (UASO 12-m) observations 
made in the BIM A Survey of Nearby Galaxies (SONG, 
iHelfer et al.ll200,^ . The combined observations resulted 
in a synthesized beam of 3.5" which projects to a linear 
scale of 75 pc. Details for the individual observations ap- 
pear in Tabled We calibrated the visibility data using 
the same techniques as employed in SONG. The visibility 
data were inverted using natural weighting to maximize 
signal-to-noise ratios. The high resolution CO observa- 
tions had slightly different pointing centers with respect 
to the low resolution and zero-spacing data, requiring in- 
version to the image domain of the high and low resolu- 
tion data separately. We then combi ned these CO maps 
using the linear combination method l|Stanimirovic et al.l 
[1999; Heifer et al. 2003). Ima ges were then cleaned us- 
ing a hybrid Hogbom / Clark / SDI algorithm with the 
residuals of the cleaned components added back in to 
avoid any flux loss. To check for possible contamination 
by continuum emission from the AGN, we imaged the 
LSB data from ^^CO and HCN. The lack of sources in 
the resulting maps places a 3(7 point source sensitivity of 
3.4 mJy at 85 GHz for objects in the galaxy and a 3.2 
mJy limit at 106 GHz. 

Although M64 has been partially mapped with the 
IRAM 30-m telescope (Casoli & Gerin 1993), the data 
are no longer available in an electronic format (F. Ca- 
soli, private communication). To provide a fully-sampled 
-"^•^CO map, we observed the galaxy with the FCRAO 
14-m telescope. We observed 7 points arranged in a 
hexagonal pattern in order to obtain a Nyquist-sampled 
map and recover all spatial frequencies. Observations of 
SiO maser stars through the night were used to check 
the pointing solution of the telescope. We observed the 
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Fig. 1. — Example of resulting spectrum from stitching pro- 
cedure. This spectrum is taken at a2000 = 12'^56™42!7, ^2000 = 
-|-21°41'18'.'5 spectra are merged around the LSR velocity of the 
galaxy, shown with the dotted line at 406 km s~^. This shows the 
general success of producing a combined high-resolution spectrum 
from different velocity configurations. 



galaxy by alternately pointing at the galaxy with two 
different pixels in the SEQUOIA 32-pixel receiver. We 
used the time that each pixel was pointed at blank sky 
as the reference position for that pixel. Vane calibra- 
tion was used to set the flux scale. We used two back 
ends in parallel: the Quabbin Extragalactic Filterbank 
(QEF), which provided 320 MHz of spectral coverage in 
64 5-MHz channels (864 km s^^ divided into 13.5 km 
s"^ channels), and the Focal Plane Array Autocorrela- 
tion Spectrometer (FAAS), which provided 80 MHz of 
spectral coverage with 0.3125 MHz channels (216 km 
s~^ divided into 0.84 km s~^ channels). Unfortunately, 
the velocity range of 216 km s~^ provided by the FAAS 
correlator is insufficient to cover the entire galaxy and 
the velocity resolution of the QEF was too coarse to 
merge with the interferometer data. Therefore, we had 
to observe the galaxy with two separate velocity config- 
urations which measured the approaching and receding 
halves of the galaxy separately. We used the QEF data 
to establish the common scaling between the two veloc- 
ity configurations. Since the FAAS and QEF follow the 
same amplifier chain, smoothing FAAS data to the QEF 
resolution should result in the same spectrum. To stitch 
together the two sets of FAAS data, we subtracted a 
linear baseline from the QEF data and used the result- 
ing spectrum as a model of signal emission in the FAAS 
data. By accounting for the expected signal using the 
QEF data, a linear baseline could be removed from each 
of the FAAS spectra separately and spectra from the two 
configurations could be merged. An example of this ap- 
pears in Figure^ Typically, the spectra had a noise level 
of 12 mK in a 4.25 km s~^ channel after correcting for a 
main beam efficiency of 0.5 at 110 GHz. As a check, we 
smoothed the final FAAS spectrum to the resolution of 
the QEF data and found there were only small variations 
resulting from higher order baselines in the data sets. 

We merged the zero-spacing information with the ^^CO 
interfe rometer data using the Linear Combination tech- 
nique l|Stanimirovic et al.llT999t IHelfer et al1l200^ . The 
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TABLE 1 

Properties of BIMA Array Observations 



Tracer 




Array 


Length 
[hours] 


[K] 


Beam Size 
["] 


^ rms 

[Jy bm"-'^ km s"-"^] 


12CO(l ^ 


0) 


B 


11 h 20 m 


610 


4.6 X 3.6 


0.24 






Cb 


10 h 57 m 


417 






13C0(1 ^ 


0) 


B 


22 h 35 m 


264 


3.3 X 2.6 


0.11 


C 


12 h 05 m 


430 










D 


3 h 35 m 


305 






HCN(1 


0) 


B 


16 h 01 m 


198 


4.0 X 3.4 


0.037 






C 


7 h 47 m 


221 










D 


3 h 58 m 


194 







^Using a 4.25 km s ^ channel. 

array data in ^^CO(l 0) are from the BIMA Survey of Nearby Galaxies 
IRee-an et idSmnt) . 



merged dataset was deconvolved using a CLEAN algo- 
rithm which terminated at the 1. barms level of the data. 
The combined dataset contains 94 ±7% of the single dish 
flux. Finally, we regridded the CO and HCN images to 
match the spatial and velocity scale of the ^■^CO data to 
facilitate pixel-to-pixel comparison, resulting in a final 
channel width of 4.25 km s^^ in all maps. 

2.2. Recombination Line Images 

To obtain recombination line emission images, we 
downloaded HST images of the galaxy in Paa (F187N), 
H-hand continuum (F160W), Ha (F656N) and ^-band 
continuum (F547M) from the HST data a rchive. The 
Pag image was presented originally by IBoker et al.l 
()1999(1 . We followed the reduction steps presented in 
that work including removal of the 'pedestal effect' in 
the NIC3 camera. However, we adopted a different 
optimal scaling between the F160W continuum image 
(Wij) and the F187N line image (Nij) for continuum 
subtraction. We minimize a goodness-of-fit parameter 
— li^ij ~ ^ith respect to the 

scaling m, where aij is the noise associated with each 
pixel and the power p is determined by 

_ r 0.5, Nij > mW,j 
\ 6, N,, <mW,j. 

This construction preferentially fits those points in the 
narrow band image that are free from Paa emission. For 
the NICMOS data, m = 0.0418. This method imphcitly 
assumes that the color of the continuum does not vary 
across the map. 

The Ha images have not been published, and were ob- 
served under HST GO Proposal #8591 (Richstone, D.). 
We used a similar procedure to that of the Paa analysis 
to extract an image containing only Ha line flux from 
the archival data. The principal differences were that we 
used median filtering among an array of images to elimi- 
nate cosmic rays. We used a F547M image of the galaxy 
for the continuum image and the F656N image contain- 
ing the line flux. The optimal scaling between these two 
images is, on average, m — 0.304 though data from each 
chip of the WFPC2 were treated separately. 

3. ANALYSIS OF THE MOLECULAR LINE MAPS 

In Figures 121 to 21 we present velocity-integrated maps 
of the molecular disk in M64. The data have been masked 



to highlight the structure of the emission in the three 
maps. Since all three data sets have the same coordinate 
axes, a single mask is generated common to all three 
data sets. An element is included in the mask if it has 
detectable emission in CO. For the CO data, the masking 
process begins by searching for all positions with emission 
(J) greater than 4 times the rms value of the noise (arms) 
in two adjacent velocity channels. These regions are then 
expanded to include all pixels with / > 2arms that are 
connected in position or velocity space by > 2arms de- 
tections to the j oint Aarm.n detection "cor e" following the 
method used in lEnerarefiola et aD l)2fl03r i. For compari- 
son, we have plotted an integrated intensity image us- 
ing a 2a„ns clip in the right-hand panel of Figure (21 
The mask generated from the CO is then applied to 
the ^'^CO and HCN emission. Using a common mask 
among the different tracers facilitates direct comparison 
of the morphology of the molecular emission. We use this 
process since emission from rare tracers will be located 
in the regions of the datacube where there is emission 
from the more common tracers. Because of this assump- 
tion, the masked maps have complex statistical proper- 
ties since the CO map contains only positive emission 
while the ^^CO and HCN emission include noise in the 
data. Nonetheless, the masking process highlights faint 
HCN emission in the disk of the galaxy which is difficult 
to detect by analyzing the HCN emission alone (Figure 
01). These faint features align with the strong CO fea- 
tures seen in Figures (21 and (31 We used Monte Carlo 
simulations to establish point source completeness limits 
with this masking algorithm for emission with a velocity 
FWHM of 2 channels. The quoted limits in Table (21 are 
the threshold for a point source to be recovered in 90% of 
the Monte Carlo trials. We use this masking process to 
examine the structure of the emission and for the rota- 
tion curve but use simpler methods for line ratio studies 
and cloud decomposition. 

The three maps share several common features. The 
emission is strongly peaked at the center of the galaxy 
with a surrounding disk. The NE half of the disk shows 
stronger emission than the SW half. The molecular emis- 
sion in the NE half of the galaxy aligns with the promi- 
nent dust lane that earns the galaxy's nickname of the 
'Evil Eye;' however, these observations show that there is 
significant molecular gas in the SW of the galaxy as well. 
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TABLE 2 

Flux Recovery for BIMA Observations 



Tracer 


Flux 


SD Flux 


Frac. 


Completeness Lmt.'^ 




Jy km 


Jy km s^^ 


Recovered 


Jy km s^-*^ 


CO 


1678 ± 23 


1845 ± 217=' 


91 ± 12% 


1.06 


i^CO 


160 ±6 


230 ± 10 


70 ± 6% 


0.35 


HON 


35 ±3 


36.5 ±5.6'' 


96 ± 18% 


0.20 



IHelfer et al.l <2003) 
^ IHelfer fc Blitl 11997bl) 

'^For recovery of 90% of point sources with velocity FWHM of 8.5 km 
s~^ (two channels). 




Fig. 2. — Masked integrated emission map of M64 in ^^CO(l — > 0) for two different masking methods. The masking method used in 
the paper (left panel) is compared to a mask generated by clipping the emission above 2arms (right panel). Comparing the two methods 
highlights the ability for the adopted method to recover low surface brightness emission. In the left panel, contours are 2^ arms for A: > 
and in the right panel, contours are 2^ cFrms for J > 3. . The levels of the contours relative to the grayscale are indicated in the colorbar 
by dotted lines. 
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Fig. 3. — Masked integrated emission map of M64 in ^^CO(l — » 
0). The contours are —2arms, —^crrms (dotted) and 2^ arms for 
fc > (solid). The levels of the contours relative to the grayscale 
are indicated in the colorbar by dotted lines. The emission appears 
significantly more clumpy than the ^^CO emission despite similar 
UV coverage. 



The disk appears to cut off sharply at a galactocentric 
radius of Rgai = 1000 pc. To quantify this, we plot the 
azimuthally averaged surface brightness distributions in 
FigureEl The azimuthal average shows a significant cut- 
off at Rgai — 925 ± 25 pc in the CO emission, well above 
the 3<Trms noise level of the observations. A similar cutoff 
is also seen in the ^^CO emission and likely exists in the 
HCN emission but is below our significance threshold. 

3.1. The Rotation Curve of M64 

M64 is peculiar because the outer H I disk {Rgai > 2 
kpc) is counterrotating with respect to the inner H I 
and stellar disk s. Studies by Braun et al. ( 1994) and 
iRix et aP l)1995j) suggest that the outer disk is result of 
a retrograde merger with a small galaxy. They propose 
that interaction between the counterrotating gas disk and 
the prograde stellar+gas disk dissipates significant angu- 
lar momentum, fueling infall to the central disk which 
accounts for the high values of the surface mass density 

of neutral gas (Rix et al 1995; Rubin, .19941 . However, 

the gas in the central region of the galaxy {Rgai ^ 1 kpc) 
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Fig. 4. — Masked integrated emission map of M64 in HCN(l — 0). 
The contours are —2arms, —io-rms (dotted) and 2*(7rms for A; > 
(solid). The levels of the contours relative to the grayscale are in- 
dicated in the colorbar by dotted lines. The emission is strongly 
peaked in the nuclear region of the galaxy, but the masking pro- 
cedure highlights faint emission, particularly in the NE half of the 
disk. 
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Fig. 5. — Azimuthally averaged surface brightness for molecular 
tracers. The averages were taken in 50 pc bins. The error bars 
represent la errors and the lines indicating upper-limits represent 
the profile that would be observed if the blank pixels contained 
emission at the Sarms level, a very conservative upper limit. The 
averages have been scaled by cos i = cos 56° to reflect face on 
values. The CO emission shows a marked truncation at Rgal = 
925 ± 25 pc, indicated by the vertical dotted line. 



appears to be in a normal, rotating disk which is decou- 
pled from the peculiar dynamics of the outer galaxy. The 
outer edge of the molecular disk occurs at the interface 
between the well-organized inner disk and the infalling 
ionized gas in the in region 1000 pc < Rgai < 1500 pc. 

Since these are the first millimeter observations to 
cover the entirety of the molecular gas disk, we exam- 
ined the dynamics of the inner region of the galaxy. 
We generated velocity fields from the masked emission 
for each of the three tracers. We analyzed the velocity 
fields using the tilted-ring algorithm first presented by 
IBegemaiil fl989) to generate a rotation curve a s imple- 
ment ed in the NEMO stellar dynamics toolkit ijTeubenI 
11995(1 . We used literature values for initial estimates of 
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Fig. 6. — Rotation curve of M64. The data agree well until the 
small numbers of points that show ^^CO emission lead to poor fits 
in the outer section of the galaxy. Typical errors in the derived 
rotation velocities vary between 0.5 and 10 km s~^ and Icr error 
bars are plotted for the ^^CO points. We also plot the radial 
velocity derived from an axisymmetric inflow/outflow model. 



the circular rotation speed, position angle, inclination, 
systemic velocity and dynamical center (B94, NUGA). 
We then iteratively fit for these parameters until a con- 
vergent mod el was determined following the procedure 
described in iWong et alJ l|2004(l . The derived parame- 
ters of the disk do not vary significantly for the molec- 
ular emission within 40" of the galactic center, so we 
forced the inclination, position, systemic velocity, and 
position angle of the rings to be equal in order to im- 
prove the estimate of the dynamical properties. The dy- 
namical center of the galaxy is at 02000 — 12*^ 56™ 43^6 
and ^2000 — 21° 40' 57'.'7 with a systemic velocity of 
Vsys = 411 ± 1 km s^^. The disk of the galaxy has a po- 
sition angle of —67° ± 3° and an inclination i = 59° ± 2° 
averaged over Rgai < 40" using 4" bins including correc- 
tions in the quoted errors for beam smearing and over- 
sampling. These dynamical properties agree well with 
those adopted by other authors. The peak of the molec- 
ular emission is slightly offset (1") from the dynamical 
center of the galax y, as is found in many disk galaxies 
l|He]fer et a]Jl200,'^ . 

The rotation curve of the galaxy is shown in Figure |3 
for each of the three tracers and the H I (from B94) . The 
rotation curves agree well for the three tracers in the in- 
ner portions of the galaxy. At large radii, the H I rotation 
curve differs from that of the molecular gas. This may 
be a physical effect arising from infalling material inter- 
acting with the molecular disk. Radial motions are seen 
in the atomic gas at larger radii {Rgai ^ 80", 1.6 kpc, 
B94), so we fit the amplitude of an axisymmetric radial 
gas flow from the velocity field of the galaxy, which is 
plotted as Vradiai in Figure El We find that there are sig- 
nificant radial motions of molecular gas for Rgai > 40" 
(800 pc), with amplitudes comparable to the motions 
found at the inner edge of the counterrotating H I disk 
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-70 km B94). In the region between 
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40" < Rgai < 60", we may be observing the interface 
between infalling gas and the rotating molecular disk. 
However, the molecular disk appears fairly asymmetric 
beyond Rgai = 40" (800 pc) and radial motions may be 
the result of elliptical gas orbits. We also confirm the 
result of M02 who noted that that the CO emission at 
the NE end of the minor axis is significantly displaced in 
velocity from the circular rotation of the galaxy. Finally, 
we see the same evidence for streaming motions in both 
CO and ^'^CO emission that is discussed in NUGA in the 
inner parts of the molecular disk. 

3.2. Line Ratios 

The CO map (Figure [^J appears to be significantly 
smoother than that of the ^^CO map (Figure Ini- 
tially, it is unclear whether this reflects differences in the 
actual emission distribution or is due to noise because 
the dynamic range in the CO map is twice as large as 
that in the ^^CO map. However, as we discuss below, 
the difference in the morphologies is real and arises from 
variations in the fraction of molecular gas found in a dif- 
fuse, translucent phasc^. 

Because the abundance of ^^CO is lo wer than that of 
CO typically by a factor of 40 to 70 l|Wilson fc RoodI 
11994(1 ■ the optical depth of the ^'^CO line is significantly 
smaller than that of the CO line. In a simple radia- 
tive transfer model through a slab of chemically uni- 
form, isothermal gas, the line ratio i?i3 = Wuco /I^iaco 
should fluctuate between the abundance ratio of the two 
species in the optically thin limit to unity in the opti- 
cally thick limit. Averaged over a po rtion of the Milk y 
Way disk, R^ = 6.7 (33.5° <b< 35°. lPo1k et ^[1981) . 
However, the value of i?i3 varies depending on the struc- 
tures observed. G iant Molecular Clouds have R13 = 4.5 
l|Polk et alJ I1988D while R13 = 10 20 in translu- 
cent, high-latitude molecular clouds l|Blitz et alJ 11984 
iMagnani et al. 1985) and R13 = 21 z fc 8 is observed in 
small molecular structures in the ISM l|KnaDD fc Bowerij 
IT9881I . At the most fundamental level, the variations 
in i?i3 reflect variations in the relative opacities of the 
two line transitions. Unfortunately, radiative transfer in 
molecular gas remains a daunting problem with solutions 
existing only for relatively simple problems. Nonetheless, 
since bound, opaque clouds show a systematically lower 
value of i?i3 than do diffuse, translucent clouds, we view 
variations in i?i3 as empirical evidence for a changing 
fraction of material in bound clouds. 

In addition to having a higher value of -R13, translu- 
cent molecular structures also exhibit a lower CO-to-H2 
conversion factor. We express the CO-to-H2 conversion 
as 



A^(H2: 

where X2 



= 2 X 10^° cm-2 • X2 



1^1200 
, 1 K km s" 



Strong fc Mattox 
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1 fo r mas s ive clouds in the Milky Way 
iDame et all |200H) and M33 
In translucent molecular 



clouds, X2 ~ 0.25 l|de Vries et~aLlll987t iMagnani et all 

^ Throughout the paper, "diffuse" gas refers to molecular gas 
that is not self-gravitating and is used in contrast with "bound" or 
"self-gravitating." Diffuse gas is often also "translucent" meaning 
Av < 3 as opposed to bound gas which is often "opaque" meaning 



l200l . i.e., CO emission from translucent structures is 
overluminous relative to the molecular column in translu- 
cent clouds when compared to more massive structures. 
Since these regions also exhibit high values of i?i3, the 
measurements of i?i3 can, therefore, be used to highlight 
regions where the assumption of X2 = 1 is questionable. 

Variations in i?i3 have been observed in e xtragalac- 
tic molecular clouds. iWilson fc Walkeil l|1994|) used the 
Owens Valley interferometer to show that i?i3 = 7.5 in a 
CMC in M33; but when the CMC and the surrounding 
gas was observed using a single dish telescope, they found 
^13 = 10. They attributed the difference between the 
i?i3 values on different scales to the presence of diffuse 
molecular gas surroundi ng the CMC wi t h i?i3 > 13.5. 
Beyond the Local Group. iPaglione et al.l l|2001fl find sig- 
nificant variations in i?i3 over the f ace of galaxi e s with 
increases in the centers of galaxies. lAalto et"al] lll995D 
find that starbursting or interacting galaxies have sys- 
tematically higher values of i?i3. These regions likely 
contain a larger fraction of diffuse molecular gas than is 
seen in the Milky Way disk. 

We use the variation of i?i3 across the face of M64 as 
proxy for whether the emission appears to rise predomi- 
nantly from highly opaque objects (i?i3 4) or translu- 
cent objects (i?i3 > 10). Since the masking process 
likely introduces systematic effects into the generation 
of moment maps, we use the unmasked data. In order 
to increase the significance of a measurement, we aver- 
aged together several spectra from annuli in the plane 
of the galaxy. We used the rotation curve from t l3.1l to 
shift each spectrum by the projected rotation velocity at 
each position which places the emission from each spec- 
trum at a common velocity. The magnitude of the shift 
is —V{Rgai) cosO smi, where 9 is the polar angle from 
the kinematic major axis, measured in the plane of the 
galaxy, and i is the inclination. This correction does not 
account for radial motions in the galaxy which are signif- 
icant at Rgai > 800 pc. This averaging process is so ef- 
fective that it permits the measurement of HCN emission 
to large galactic radii, even though little emission can be 
clearly identified in the data cubes. Examples of the av- 
eraging appear in Figure|7|for two annuli: Rgai < 100 pc 
and 500 pc < Rgai < 600 pc, which show substantially 
different line ratios, particularly in HCN. The spectrum 
for Rgai < 100 shows evidence for the HCN/CO line ratio 
varying over the line with the wings of the line showing 
comparatively more HCN emission than the core of the 
line. The wings of the line are generated from region at 
Rgai ~ pc where the velocity gradient across the syn- 
thesized beam is the largest. The larger line ratio in the 
wings suggests that the molecular gas is densest in the 
region near the AGN (c.f., NUGA). 

We measured the line ratio as a function of position 
across the galaxy. We integrated the emission over the 
derived spectra to measure the flux in annuli of galac- 
tic radii 100 pc wide. We then took ratios of the fluxes 
and plot the results in Figure |H1 The value of i?i3 
changes significantly over the face of the galaxy with 
most of the molecular disk having a value of R13 ~ 6 
comparable to the disk of the Milky Way. The nuclear 
region {Rgai < 200 pc) and beyond Rgai — 800 pc show 
an increased value of -R13. The material at small galac- 
tocentric radius may be affected by the presence of low- 
luminosity AGN and a nuclear starburst l|Ho et alJll997D 
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Fig. 7. — Two examples of the spectra derived from the averaging 
process used to measure Une ratios. The top panels shows the 
emission averaged over Rgai < 100 pc and the bottom panel shows 
the average for 500 pc < Rgai < 600 pc. Emission from ^^CO has 
been scaled up by a factor of 6.5 and emission from HCN has been 
scaled up by a factor of 12. The line ratios change significantly 
over the face of the galaxy. 



which results in an increase in the kinetic temperature of 
the molecular gas. Such an increase in the temperature is 
thought to produce a n increased value of J? i .■^ near the nu- 
clei of many galaxies iPael ione et al.l l20011 and is consis- 
tent with the enhanced C0( J = 2 ^ 1)/C0(J =1^0) 
ratio observed at Rgai < 100 pc in the NUGA study. At 
large radii, the material may be disturbed by infalling 
material from the counter-rotating atomic disk beyond 
Rgai = 1000 pc because the increase in Ri^ is seen in gas 
where radial motions are observed f H3.1|l . The change in 
i?i3 is probably because the fraction of gas in a diffuse, 
translucent phase is higher in this region than in the 
remainder of the galaxy. In both of these regions, the 
standard CO-to-H2 conversion factor is suspect, since it 
is derived from gas where most of the molecular emission 
arises from cold, bound, opaque GMCs. 

Near Rgai ~ 0, the ratio Whcn/W^co is comparable to 
that se en in the nucleus of the Milky Way ijHelfer fc Blit j 
Il997aj) On average, the value of W^hcn/Wco over the en- 
tire galaxy is comparable to the value seen in the high 
column density regions of GMCs in the Solar neighbor- 
hood or the inner disk of the Milky Way. The relatively 
high values of this line ratio imply that a substantially 
larger fraction of the molecular gas is found in high den- 
sity regions than is typical for the Solar neighborhood 
where Wkcn/Wco < 0.01 ifHelfer W\tM^c,7^ . 



3.3. Fourier Analysis of the CO and CO Line Ratios 

The integrated intensity images of the CO and ^^CO 
emission look qualitatively different; the ^^CO emission 
looks significantly "clumpier." Since the "'^'^CO emission 
should be a more faithful tracer of molecular column den- 
sity than CO emission, the clumps seen in Figure O sug- 
gest that there are clouds present in the data. However, 
the ^^CO map also has lower dynamic range and the ef- 
fects of masking may produce the difference in structure. 
In order to determine whether the ^'^CO map is actually 
clumpier, we measured the variation in line ratio as a 
function of spatial scale using Fourier techniques. Since 
bound clouds have a fixed spatial scale {£ ^ 100 pc in 
the Milky Way), emission arising from them should have 
most of its power around this scale. In contrast, emis- 
sion from diffuse clouds is not limited to a fixed scale and 
should contribute at all spatial scales. Since the line ra- 
tio i?i3 is different in diffuse vs. self-gravitating clouds, 
the ratio of power in the CO map to that of the ^"^CO 
map should vary as a function of spatial scale. 

We generated an integrated intensity map of the galaxy 
by shifting the emission in the spectra to a common ve- 
locity in order to eliminate the effects of galactic rotation. 
We integrated over a window 106.25 km s~^ wide (25 
channels) centered on the emission line. This generated 
integrated intensity maps with improved signal-to-noise 
without introducing any erroneous spatial filtering from 
masking the data. The resulting maps contained > 80% 
of the emission found in the original maps with the miss- 
ing emission in the nuclear region where the linewidth 
exceeds the width of the window. Next, we Fourier trans- 
formed the integrated intensity maps and averaged the 
power in elliptical annuli so that each annulus represents 
a range of scales in the deprojected image. We then sub- 
tracted the noise power from the data using the trans- 
form of noise images generated from emission-free regions 
of the data cube. The resulting power as a function of 
scale should represent only the power from line emission. 
In Figureiniwe plot the results of the analysis, scaling the 
^■^CO image up by the ratio of the fluxes at small scales 
(4.1) to facilitate comparison. The power spectra have 
the same shape at large scales {£ > 200 pc) and small 
scales (75 < £ < 200 pc) but they have different ratios 
of powers in these two regimes: 6.5 ± 0.3 for £ > 200 pc 
and 4.1 ± 0.1 for 70 < £ < 200 pc. The value of i?i3 
at small scales is comparable to that seen in Galactic 
GMCs, and emission on these scales is from structures 
that may be the analogues of Galactic GMCs. Compared 
to the CO map, the ^■^CO map has less emission at large 
scales where only diffuse gas should contribute. Hence, 
^•^CO emission from clouds will be less confused by the 
presence of diffuse gas compared to CO emission. 

The flux ratios derived in this analysis do not agree 
with the ratios calculated for the galaxy as a whole (Fig- 
ure for two reasons: (1) much of the emission with a 
large value of i?i3 is located in the center of the galaxy 
where the velocity range is smaller than the linewidth so 
some high-i?i3 emission is not included; and (2) incom- 
plete subtraction of the noise power will skew the ratios 
to lower values. This is because the noise power, relative 
to the signal, is larger in the ^'^CO emission. If there 
is incomplete subtraction of the noise, it will contribute 
more to the power in the "'^^CO emission than to the CO 
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Fig. 8. — Line ratios between molecular gas tracers as a function of galactocentric radius. Over most of the disk, the CO to ^^CO ratio 
has a values typical of the Milky Way disk (~ 6) but shows a significant increase at small and large values of Rgai- For Rgal > 800 pc, 
there are significant radial motions which may introduce systematic errors into the measurement technique. These increases may indicate 
significant amounts of difi'use molecular gas. The HCN/CO ratio is larger than seen in the GMCs in the Solar neighborhood, implying the 
average volume density of molecular gas is higher in M64 than it is in the Milky Way. 
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Fig. 9. — Power spectrum of CO and ^'^CO images as a function 
of spatial scale. The figure compares the power spectrum of the 
original CO data and the ^^CO image after being scaled by ratio of 
the fiuxes at small spatial scales (4.1). The distribution of emission 
agrees well on the small scales but the line ratio at large spatial 
scale scales is larger (6.5). The dilference in line ratios as a function 
of scale implies that the clumpy structure in Figure |21 is real and 
the line ratios are comparable to those seen in GMCs in the Milky 
Way. 



emission. 

The observed difference in flux ratios as a function of 
scale persists even when the dynamic range of the CO 
image is reduced to match that of the ^'^CO image. We 
conclude that the different morphologies in Figures[21and 



El represent a real difference in the structure of the emis- 
sion from the galaxy and are not the result of the mask- 
ing process. Furthermore, it appears that there is sig- 
nificant contribution to the molecular emission from a 
diffuse, translucent component of the molecular ISM. By 
analogy with molecular emission in the Milky Way, it is 
likely that it has a conversion factor X2 < 1. Despite the 
lower signal-to-noise, the "'^'^CO map should be a better 
tracer of the molecular emission from bound clouds in 
M64. We return to the question of diffuse gas in US. II 

4. DECOMPOSITION OF THE MOLECULAR EMISSION 

To answer the fundamental question of this study, 
namely 'Are there GMCs in M64?', we decompose the 
blended molecular emission into clouds and then com- 
pare the properties of these clouds with Local Group 
GMCs. For this purpose, we focus on the ^'^CO emis- 
sion since this tracer will be less contaminated by emis- 
sion from diffuse molecular gas. There are three prin- 
cipal methods for isolating structures in blended molec- 
ular emission (1) partitioning of the cube by eye (e.g., 
IWilson fc ScovnafToOfi IWright et aLllTool . (2) "clean- 
ing" the data set by iteratively fitting and subtracting 
three dimensional Gau ssians pioneered with the GAUSS- 
CLUMPS algorithm bv |Stutzki fc GiTstCTl EMh) (3) par- 
titioning of the emission by assigning it to the "closest" 
local maximum first developed in the CLUMPFIND al- 
gorithm by WiUiams ct al. (1994'). We adopt the latter 
approach for M64 since it docs not require known func- 
tional forms for the ^^CO distribution of the clouds and it 
works directly with the emission present in the datacube. 
The original CLUMPFIND algorithm was implemented 
for the high signal-to-noise regime on single dish data. 
We have redevelop ed the algorithm using the methods of 
iBrunt et al.l l)2003j) to make the method more robust in 
the presence of noise. In addition, we have modified the 
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distance metric in the algorithm to be apphcable to in- 
terferometer data which is oversampled in the spatial di- 
mension. A detailed discussion of the algorithm appears 
in Appendix including benchmarking information for 
the choice of parameters. The algorithm can be applied 
to noisy, interferometric data with minimal variation in 
the properties of the derived clouds as a result of small 
changes in the free parameters of the algorithm. 

Despite the modification of the algorithm to minimize 
the effects of of noise, decomposition of the ^'^CO masked 
data used in Figure 3 was hampered by including low 
significance emission and noise. To improve the sensi- 
tivity of the datasets, we smoothed all three datacubes 
to a channel width of 8.5 km s~^. We regenerated a 
mask for the CO data using the same method described 
in namely we required / > 2arms in a pixel that 
is connected to / > -iarms detections in two adjacent 
channels. However, for the ^'^CO and HCN data, we re- 
quired / > 2arms in two adjacent channels where there 
was detectable CO emission. The threshold was cho- 
sen so that < 1 false detection would be included in the 
^■^CO and HCN masks. While this method discriminates 
against low significance emission, it does a far better job 
of selecting only real emission in the weaker tracers of 
molecular gas. 

We applied the decomposition algorithm to the molec- 
ular emission in all three tracers. We also applied the 
algorithm to the CO data at the full velocity resolution 
because of the high signal-to-noise of the original CO 
data. In Figure E| we display the masked, integrated 
intensity images for ^'^CO and HCN, including the lo- 
cations of the clouds into which the emission is decom- 
posed. The algorithm isolates 87 clouds in the ^^CO 
emission and 42 clouds in the HCN. Only a fraction (25 
and 9, respectively) of these clouds are resolved in both 
position and velocity space and these clouds are indi- 
cated with stars. The decomposition of the CO data at 
full velocity resolution yields 44 clouds. The algorithm 
identifies several clouds that appear in all three tracers 
with similar properties. 

4.1. Cloud Properties 

We measure the macroscopic properties of these clouds 
using intensity- weighted moments of the emission distri- 
buti on, following the method s of iSolomon et a 
and iRosolowskv et alJ l|2003j) to facilitate comparison 
with Local Group GMCs. Using these methods, we de- 
termine the mean position and line-of-sight velocity of 
each cloud. In addition, we measure three macroscopic 
properties for the clouds: the velocity FWHM (AV^), 
luminosity (L) and deconvolved effective radius (Re)- 
We have corrected the derived values for these proper- 
ties to account for the fact that the data are clipped at 
TcHp = 2arms] and therefore, the properties are not mea- 
sured over the same dynamic range as for Local Group 
clouds. To facilitate comparison with ot her data sets, 
we ad opt th e Gaussian correc tion used in lBolatto et al.l 
l|2003D and lOka et alJ l|2001l^ . The derivation of the 
macroscopic cloud properties is given in detail in Ap- 
pendix IH Table El is a listing of the cloud properties for 
the "CO clouds. 

4.2. Comparison with Milky Way CMCs 



In the Local Group, the macroscopic properties of 
molecular clo uds are r elated by power law relationships, 
first noted by 'Larson' ('1981'). In particular, the study of 
.Solomon et al. ( 1987) found 



AV 

km s^ 

il2CO 



= 1.7(:^ 

pc 



0.5 



and 



K km s ^ pc ^ 



= 1.8 



AV 

km s^ 



We use these values for comparison with the results of 
this study since we have attempted to replicate their 
measurement methods and systematic effects. To com- 
pare the derived ^^CO luminosities to the luminosities 
of S87, which are measured in ^^CO, we scale their re- 
lationship down by a factor of i?i3 = 6.7, the global 
CO/^'^CO line r atio, in the inner disk of the Milky Way 
ijPolk et alJl988il . In Figure^2 we plot the relationships 
between these macroscopic properties for the clouds in 
M64. We have excluded clouds with linewidths larger 
than AV > 60 km s~^. These clouds contain multiple 
peaks and appear to be the blend of two smaller clouds. 
These clouds are found near the center of the galaxy 
where the large velocity gradient makes decomposition 
of the emission difficult. 

Several features are immediately apparent from Fig- 
ure ^2 First, the properties of the clouds are signifi- 
cantly different from the relationships derived by S87. 
Second, the clouds have large luminosities compared to 
the Local Group clouds. In S87, very few clouds are 
massive enough to have a ^'^CO luminosity larger than 
10^ K km s~^ pc^. For reference, the mass implied by 
the CO luminosity has also been plotted in Figure ITTI 
assuming i?i3 = 6.7 and X2 = 1, which gives a ^^CO- 
to-H2 conversion factor equivalent to that derived us- 
ing dust extinction by Lada et al. (1994)- To compare 
the Milky Way and M64 clouds in a more quantitative 
fashion, we fit hnear trends between log(i?e), log(Ay) 
a nd logfLi.^no) for the -'^'^ CO clouds using the method 
of lAkritas fc Bershadvl l)1996f >. which is appropriate for a 
population with intrinsic scatter and measurements with 
errors in both parameters. We find that 

log (AV) ^ (-0.6 ± 0.4) + (1.0 ± 0.3) log (Re) and 

log (ii3Co) = (2.2 ± 0.6) + (2.2 ± 0.4) log (AV) . 

The index for the linewidth-size is marginally different 
from the value for the Local Group (1.0±0.3 vs. 0.5), and 
the luminosity-linewidth is significantly different from 
the Local Group trend (2.2 ± 0.4 vs. 5.0). These fits 
indicate that the clouds identified in the M64 data are 
not like molecular clouds seen in the disk of the Milky 
Way. However, the luminosity of these clouds implies 
that most of the clouds are significantly more massive 
than any clouds an a lyzed in the S87 study, and indeed 

argue that no clouds more 
massive than M > 6 x 10^ Mq will be found in the 
Milky Way disk. Thus, we are comparing the proper- 
ties of these clouds with an extrapolation of Local Group 
clouds to high mass. It is noteworthy, however, that the 
low luminosity clouds in M64 have comparable luminosi- 
ties, sizes and linewidths to the most massive clouds in 
the Milky Way. Thus, the M64 clouds may indicate a 
change in the scalings among CMC properties that oc- 
curs at high mass (or high average surface density). We 
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Fig. 10. — The results of the decomposition algorithm applied to the ^"^CO and HCN data sets. The locations of the clouds are overlaid 
on a grayscale integrated intensity maps of the galaxy. Locations of the clouds are indicated by stars (for clouds that are resolved in both 
position and velocity) and dots (unresolved in either position or velocity). The algorithm identifies 87 clouds in the ^^CO data, of which 
25 are resolved, and 42 clouds in the HCN data, of which 9 are resolved. The effective radii (-Re, see Appendix^} of resolved clouds are 
indicated by the circle drawn around them. Several of the ^^CO clouds are also seen in the HCN emission with similar derived properties. 



TABLE 3 

Properties of the resolved -'^'^CO clouds 



Number 


Position^ 


Vlsr 


Re 


AV 


^13CO 


pb 




("/') 


[km s--"-] 


[pc] 


[km s-^] 


[10-4 K km s-i pc2] 




1 


(29,-7) 


266.1 


176 


28 


87.4 


4.7 


2 


(10,-4) 


292.8 


115 


48 


79.7 


5.1 


3 


(-4,13) 


473.5 


174 


68^= 


67.0 


3.6 


4 


(-16,7) 


569.3 


146 


26 


56.6 


4.1 


5 


(-64) 


539.1 


89 


39 


55.1 


4.2 


6 


(-1,-4) 


401.6 


136 


73^^ 


54.8 


3.8 


7 


(-5,0) 


491.8 


104 


46 


53.4 


3.6 


8 


(10,15) 


399.0 


144 


38 


42.0 


4.0 


9 


(28,1) 


298.0 


68 


21 


31.8 


4.5 


10 


(2,0) 


345.1 


50 


35 


23.4 


3.5 


11 


(19,8) 


349.8 


81 


23 


21.5 


3.2 


12 


(22,3) 


318.7 


78 


19 


16.3 


3.4 


13 


(-45,7) 


550.6 


86 


18 


15.8 


3.3 


14 


(-19,24) 


506.4 


99 


25 


13.4 


2.2 


15 


(-37,27) 


544.0 


53 


18 


13.4 


3.8 


16 


(-16,1) 


541.8 


64 


28 


12.9 


3.1 


17 


(-36,-2) 


519.6 


95 


16 


10.7 


1.8 


18 


(-50,20) 


581.3 


115 


16 


10.3 


2.3 


19 


(-29,27) 


520.6 


64 


16 


9.6 


3.5 


20 


(-5,17) 


493.3 


90 


19 


8.2 


1.7 


21 


(-9,-18) 


404.6 


50 


10 


7.2 


2.5 


22 


(-31,18) 


546.5 


48 


20 


7.0 


2.2 


23 


(0,-21) 


380.3 


62 


23 


5.4 


2.4 


24 


(-27,-2) 


509.4 


51 


12 


4.9 


1.9 


25 


(30,-21) 


264.1 


43 


16 


4.2 


1.7 



''Given in arcseconds relative to the center of the galaxy at 
02000 = 12^ 56™ 43?6 and ^2000 = 21° 40' 57'.'7 

''The ratio of the peak antenna temperature to the clipping level for the clouds, 
which determines the magnitude of the Gaussian correction applied (see Ap- 
pendix |bJ. 

'^This cloud is a blend of two peaks of emission which cannot be separated 
objectively by the algorithm. 



GMCs in M64 11 




Fig. 11. — Plot of the linewidth— size and luminosity-linewidth relationship for clouds in M64. The left-hand panel plots the relationship 
between cloud sizes and lincwidths from the analysis of all three tracers. The right-hand panel plots the luminosity as a function of 
linewidth for the CO and ^^CO clouds. The luminosities of the CO clouds have been scaled down by 6.7 to compare them with the ^"'CO 
clouds. The derived properties of the clouds are different from the corresponding Local Group trends plotted as dashed, gray lines. The 
M64 clouds appear to have smaller sizes or larger linewidths (or both) than clouds in the Local Group. For reference, the right-hand axis 
of the luminosity-linewidth plot indicates the luminous masses of the clouds assuming X2 = 1 and R13 = 6.7. 
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cannot determine whether the apparent differences with 
respect to the S87 trends are an environmental effect or 
are intrinsic to high mass molecular clouds. 

While the analysis has focused on the ^^CO clouds, 
Figure^Jshows that the clouds generated for other trac- 
ers have similar behaviors. Scaling the luminosity of the 
CO clouds down by 6.7 to place them on a similar scale 
as the ^^CO clouds produces a luminosity-linewidth rela- 
tionship statistically indistinguishable from the relation- 
ship for the ^^CO clouds. In contrast, the sizes of the CO 
clouds appear significantly larger than the ^'^CO clouds. 
A two-sided KS test indicates that the distributions of 
the CO and ^^CO cloud radii are significantly different 
{Pks > 0.9993), while the distributions of the linewidths 
are indistinguishable. We suspect that the larger CO ra- 
dius results from the algorithm including the diffuse CO 
emission that surrounds the clouds in the definition of 
the CO clouds. To confirm this, we re-cataloged the CO 
emission clipping the map at 5arms, producing roughly 
the same dynamic range as the ^'^CO map. The high clip- 
ping level eliminates the low surface brightness, diffuse 
CO emission. Our inclusion of a Gaussian extrapolation 
to a common clipping level facilitates direct comparison 
between the resulting catalogs. The clipped catalog con- 
tained 33 clouds; these clouds have sizes comparable to 
the "'^'^CO clouds, and a two-sided KS test found the dis- 
tributions of both the sizes and the linewidths to be in- 
distinguishable between the two tracers. The presence of 
diffuse emission can thus significantly affect the sizes of 
clouds traced by CO emission alone. 

A major concern with the decomposition is whether 
the structures that the algorithm isolates really do rep- 
resent discrete physical entities or whether the "clouds" 
are really blends of several small clouds that cannot be 
decomposed with the resolution of these observations. 
We suggest that this is not the case and that we are 
cataloging physical objects for several reasons. First, we 
tested the hypothesis that the clouds are blends of several 
GMCs like those found in the Milky Way by analyzing 
simulated data sets. For each data set, we constructed 
a population of 50 molecular clouds with a total mass of 
2 X 10^ Mq and individual masses drawn randomly from 
a mass distribution with dN/dM cx M"-'^-^. For each 
cloud in the simulation, we calculated the radius and 
linewidth that would be expected based on the study of 
S87. We randomly distributed these clouds in a data 
cube so that they would appear as a cloud with —150 
pc and AV = 40 km s^^, which is appropriate for a typ- 
ical massive cloud in our data set (Figure ^J. We then 
convolved the simulated cube to the resolution of the ob- 
servations and added noise to produce a peak signal-to- 
noise ratio of 8. We processed each simulated cube with 
the same masking and decomposition algorithm used in 
the data analysis. The algorithm identifies substructure 
within the data cubes in 90% of the trials, leading us to 
conclude that the clouds are not likely to be blends of a 
CMC population like that seen in the disk of the Milky 
Way. 

The second piece of evidence that the decomposition 
products represent real objects is that the linewidths of 
the observed clouds, while large, are still significantly 
smaller than would be expected if the emission were com- 
posed of several smaller clouds orbiting in the galactic 



potential. We estimate the linewidth due to shear as 



shear 



A^'(a^o,yo), 



where I{xi,yi) is the integrated intensity at each posi- 
tion in the cloud, (v — Vr){xi,yi) is difference between 
the mean velocity and the line-of-sight projection of 
the rotation velocity at each position in the cloud, and 
AV'^{xo,yo) is the linewidth of the emission distribution 
at the cloud center. The final term accounts for broaden- 
ing due to beam-smearing and the intrinsic gas velocity 
dispersion. The measured linewidth for the clouds {AV 
in Table inil is smaller than AVshear by a factor of 2.3 on 
average and does not change significantly with radius, 
even in the outer galaxy where the shear decreases. 

Estimates of the cloud masses using the virial theorem 
agree well with estimates using the ^'^CO emission, giving 
a third piece of evidence that we are identifying clouds. 
In particular, if we adopt the virial mass estimate used 
by S87, 



/ Mvt \ 



= 187 



AV 

km s^ 



Re 

pc 



(1) 



we find good correspondence between the virial mass esti- 
mates and luminous mass estimates as indicated in Fig- 
ure El A linear fit between the log of the two mass 
estimates gives 



log(MvT) 



-1.5 ± 1.2) + (1.2 ± 0.2) log(Mco), 



consistent with a direct correspondence between the mass 
estimates over an order of magnitude. This good agree- 
ment again suggests that the clouds are discrete ob- 
jects, since the application of the virial theorem to un- 
bound groups of molecular clouds gives virial masses 
significantly in exce ss of the mass estimated from the 
CO luminositv (e.g .lRand fc KulkarillToM lRa^lT995l: 
lAllen fc Leaueuxll993|) . We have also used four alt ernate 
meth ods to measure the virial parameters ( Heislcr 'et alJ 
IT981 and find no significant differences between the var- 
ious methods. 

If the clouds represent self-gravitating objects, we can 
use the relationship between the luminous mass estimates 
and dynamical mass estimates to derive conversion fac- 
tors for the tracers. Assuming the clouds are virialized, 
we perform a least-squares fit between log(MvT) and 
log(Lco), fixing the slope at unity to determine the scal- 
ing between these two tracers. Using these fits, we derive 
CO-to-Il2 conversion factors: 



N(H2) = (1.8 ±0.2; 



X 10^^ cm"2 



W13. 



CO 



A^(H2) = (3.2±0.2)xl02° cm-2-T4^i2co or X2 = (1.6±0.1), 

where the surface brightness W is expressed in K km s^^. 
If we assume the clouds are marginally bound instead of 
being virialized, we derive conversion factors that are a 
2 times lower than these values. Given that the derived 
con version factors bracket the measured value s for both 
CO HDame et a,1.l200HlStToTig fc MattoxHaQfTT and ^^CO 
ijLada et alJll994|) . it appears that dynamical estimates 
of the mass agree well with masses estimated from CO 
luminosity. Significant external pressure or the presence 
of magnetic fields will reduce the value of A/vt, closing 
the small gap between the luminous {X2 ~ 1, R13 — 6.7) 
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Fig. 12. — Relationship between virial mass estimates and lu- 
minous mass estimates of the molecular gas for both ^^CO and 
CO clouds. The luminous mass estimates assume X2 = 1 and 
R13 = 6.7. The good correspondence between mass estimates for 
over an order of magnitude in mass implies that the clouds iso- 
lated by the algorithm are discrete physical entities, illustrated by 
a linear fit between the log(MvT) and log(Mco)- The locus of 
marginally bound objects for X2 = 1 and R13 = 6.7 is shown as 
the dot-dashed gray line. 



and virial mass estimates l|McKee Sz Zweib'eilll992f) . In 
light of this, we adopt X2 = 1 for the CO obse rvations 
and the conversion factor of L ada etld] ljl994|) for the 
^•^CO. Under the assumption of virialization, we also 
derive 
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(90± 10) 



L 



HCN 



K km s ^ pc 



implying HCN is brighter relative to CO in these clouds 
than is seen in Solar neig hborhood GMCs ( Heifer & Blit^ 
^97a). However, HCN/CO ratio is similar to gas in the 
nuclear disk of the Milky Way where the mean density 
is higher. 

With mass measurements of the molecular clouds, we 
can convert these estimates to surface densities of the 
clouds. We plot the surface density of molecular gas 
in the clouds as a function of cloud mass in Figure 1131 
On average, the surface densities of the clouds are larger 
than are seen in the Milky Way with a mean value of 
(250 ± 15) Mq pc^^, significantly highe r than the v alue 
of 100 Mp, pc^^ seen in the Milky Way ( lBlitzlll993D and 
M33 ( Roso lowskv ct al. 2003). In addition, the surface 
density is not constant with mass, which must be true if 
the clouds in the galaxies are self-gravitating but do not 
follow the same linewidth-size relationship as is observed 
in the Local Group. We find that 

f0.7±0.2 



jH2 



CX 



meaning that high mass molecular clouds are systemati- 
cally smaller and denser than would be expected from 
the extrapolation of trends seen in the Local Group. 
Moreover, an increased surface density of molecular gas 
also implies that the internal pressures of the molecular 




Fig. 13. — The scaling of molecular surface density (Eh2) with 
luminous mass Mi^co- The molecular clouds in M64 have signif- 
icantly higher surface densities than do clouds in the Milky Way. 
Moreover, the surface density is not constant with mass, as is nec- 
essary if the clouds are self-gravitating but follow a different size- 
linewidth relationship from clouds in the Local Group. 



clouds are larger than are typical in the Local Group. 
We estimate the internal pressures of these clouds: 



= 212 cm-3 K 



M 



Re 

pc 



/ AV 



\ km 1 



•(2) 



We find a mean internal pressure of 10® K cm"'^, an or- 
der of magnitude larger than the value inferred for Milky 
Way GMCs (Blitz 1993). However, an important quan- 
tity for comparison with the internal pressure is the mean 
exter nal pressure exerted on th e cloud by the ISM of 
M64. IBlitz fc Rosolowskvl lj2004r formulated an expres- 
sion for the midplane hydrostatic pressure of the disks of 
galaxies: 



= 272 cm"3 



-'gas 



km s 



Mq pc-' 
-0.5 



M, 



o pc~ 



0.5 



(3) 



where Vg is the 1-dimensional velocity dispersion of the 
gas, /i, is the stellar scale height, and is the stel- 
lar surface density. We assume a gas velocity disper- 
sion of Vg — 8 km s~^, corr esponding t o the linewidth 
of forbidden hue emission in IRix "eTaLl 1IT995I) . We de- 
rive an upper limit for the pressure by assuming and h^: 
of 30 pc, typical of disk galaxies (Blitz fc Rosolowskvl 
HO^. We measure E* deri ved from the 2MAS S Large 
Galaxy Atlas image of M64 llJarrett et alJl2003|) assum- 
ing Mk/Lk = 0.5 l)Bell fc de Jond 120011) . Finally, we 
take Eg = Eh2 + Ehi- We calculate Eh2 from the ^^CO 
surface dens ity profile and use Srt = 10 Mq pc~^ as 
measured in iBraun et alJ 1)19941) . With these data, we 
compare the internal pressures of molecular clouds to 



14 



Rosolowsky & Blitz 



the external pressure at the position of the cloud and 
find that the clouds are overpressured relative to the ex- 
ternal medium by at least a factor of 2. 

If we refine this estimate, the difference between Pint 
and Pext grows larger. Of particular note, this esti- 
mate has included the contribution of the gas in self- 
gravitating molecular clouds in the contribution to the 
external pressure. These clouds likely do not contribute 
to the hydrostatic support of the ISM. Since this gas 
dominates the local column density, the actual external 
pressure is likely significantly smaller than the above es- 
timate. If we use the estimate that 25% of the molecular 
gas is found in the diffuse component f tl5.1(l and assume 
that self-gravitating clouds do not contribute to the hy- 
drostatic support of the ISM, then the clouds have in- 
ternal pressures higher than the local ISM by a factor of 
12, similar to the margin by which Solar neighborhood 
GMCs are overpressured with respect to the local ISM 
(^litz 1992|. Recent star formation activity will decrease 
the mass-to-light ratio and therefore the mass estimate of 
the stellar disk, also reducing the external pressure. The 
central concentration of mass in GMCs tends to increase 
the internal pressures. Finally, the stellar scale height is 
likely larger than the 300 pc typical of galactic disks be- 
cause of the presence of a significant bulge in the galaxy 
across the inner ~ 300 pc (MoUcnhoff & Hcidt 200lT). 

Because the clouds are significantly overpressured with 
respect to the ambient medium, we take this as a final 
piece of evidence that the clouds are self-gravitating and 
that the luminosity traces molecular mass with conver- 
sion factors comparable to Milky Way values. Adding 
this to the agreement between the luminous and dynami- 
cal mass estimates and the small linewidths of the clouds, 
there is good evidence that molecular clouds identified by 
this method are gravitationally bound. The masses of 
many clouds are significantly larger than are seen in the 
disk of the Milky Way, wher e galactic tides set a n upper 
limit to the mass of a GMC (IStark fc BlitdlTOTl . There 
is significant shear in the inner disk of M64 which may 
disrupt the clouds so we measured the tidal acceleration 
across each of the GMCs. We find that the clouds in Ta- 
ble ^ar^^narginaJfysta^ against tides using the meth- 
ods of lStark fc Blit^ ifToTl . All clouds at Rgai > 400 pc 
are very stable against tidal disruption. We conclude 
that these clouds are the analogues of the GMCs seen 
in the Milky Way, but have significantly different macro- 
scopic properties. 

4.3. Robustness of the Results 

Several assumptions and choices have been incorpo- 
rated in the conclusion that the molecular clouds ob- 
served in M64 are self-gravitating structures, analogous 
to high mass Giant Molecular Clouds. However, the prin- 
cipal results of our study do not change even if we relax 
these assumptions. 

The largest influence on our cloud properties is the 
choice of the decomposition algorithm. The algorithm we 
developed has been designed to be robust against small 
variations in its parameters. Although we justify the 
choice of free parameters in Appendix ^ we have re- 
cataloged the ^'^CO emission while varying contouring 
levels and values of Tuniq (the free parameters of the 
algorithm). We found all changes in the cloud properties 
were small and well within the uncertainties. 



Because the clouds are comparable to the beam size, 
there is significant concern that the objects we catalog 
are the blends of several smaller clouds. We have ar- 
gued in M.2\ that this is not the case, based on simu- 
lations, pressure arguments, dynamical mass estimates 
and cloud linewidths. Several of our results rely upon 
accurate decomposition of the emission. However, the 
most important of our results remain even if the clouds 
are blends of smaller structures. Namely, the surface 
(and volume) densities of the small clouds must be even 
higher than measured for the M64 clouds, implying they 
are even more over-pressured with respect to the ambi- 
ent ISM; the estimate of the internal pressure (Equation 
121) becomes a lower limit. Moreover, much of the fol- 
lowing discussion (SjSJ is independent of how the clouds 
are decomposed. We also emphasize that the molecular 
structures cannot consist of Milky Way GMCs overlap- 
ping along the line of sight since Milky Way GMCs clouds 
are not dense enough to exist as discrete structures in the 
high-pressure molecular disk of M64. 

The derived properties have been corrected for the fi- 
nite signal-to-noise ratio in the dataset using a Gaussian 
extrapolation (see Appendix El- We have tested the 
validity of the Gaussian correction to the cloud prop- 
erties by fitting three-dimensional Gaussians to the re- 
sults of the decomposition. We find that the extrapo- 
lated values agree well with the fits. If we repeat the 
analyses of ii4.2l without the corrections to their prop- 
erties, we find (1) the indices of the linewidth-size and 
luminosity-linewidth relationships are still significantly 
different from the Local Group values, (2) the luminous 
and virial mass estimates still agree well, (3) the mean 
surface density of the clouds is (450 ±20) Mq pc~^, even 
higher than our adopted value and (4) the clouds are 
still significantly overpressured with respect to the ISM 
of M64. 

Finally, the distance to M64 is somewh at uncertain 
with a recent estimate of I? = 7.5 Mpc ijTonrv et alJ 
120011) . Adopting a larger distance will change the val- 
ues of the luminosity (as D'^) and the sizes (as D) but 
leave the linewidths unaffected. The indices of the scal- 
ings between the macroscopic properties will be invariant 
and remain different from the Local Group. The clouds 
will appear to have higher luminous masses relative to 
their virial masses {M^^/Mco oc D~^) and will still ap- 
pear bound, though the derived conversion factors will 
be smaller. The surface densities of the molecular clouds 
do not vary with the adopted distance and will be signif- 
icantly higher than Local Group values. While changing 
the adopted distance will alter the measured cloud prop- 
erties, the main conclusions are robust: the most massive 
clouds in M64 appear to be self-gravitating structures 
with macroscopic properties significantly different from 
an extrapolation from Local Group GMCs. 

5. DISCUSSION 

In this section, we relate the results of this study to 
molecular gas on galactic scales. We estimate the con- 
tribution of diffuse molecular gas to emission from M64 
and examine the dust content and star-forming proper- 
ties of the molecular clouds. Clarifying the properties of 
these high mass molecular clouds helps to link our under- 
standing of star formation on small scales in the Milky 
Way to starbursts and ultraluminous infrared galaxies 
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(ULIRGs). 

5.1. The Fraction of Diffuse Gas in M64 

Given the large CO to ^^CO line ratio {R13) observed 
in some portions of M64, we suspect that a significant 
contribution of the emission arises from diffuse, translu- 
cent molecular gas. The resolved GMCs in Table |31 
comprise ^50% of the ^^CO emission from the galaxy. 
We cannot distinguish how much of the remaining emis- 
sion arises from GMCs below our completeness limit or 
whether the emission arises from diffuse molecular gas. 

Since emission from diffuse gas has a substantially dif- 
ferent value of i?i3 from that found in the GMCs, we 
can use our observations to establish the ratio of CO 
emission coming from diffuse gas vs. bound objects 

{F=Wd^ff/WGMc)■■ 

p _ Rdiff Ravg — RgMC 
RgMC Rdiff — Ravg 

where Rdiff is the line ratio in diffuse gas, Rgmc 
is the line ratio in GMCs and Rav „ is the average 
of the line ratio over a given region liPolket alJlTosl 
iWilson fc WaIkCTllT99l . Within the 25 clouds identified 
in the ^'^CO emission, the line ratio i?i3 = 3.7, which 
is typical of the peaks of emission from GMCs in our 
own galaxy (Polk et al. 1988). This suggests that the 
physical conditions averaged over these GMCs are com- 
parable to the highest column density regions of GMCs 
in our own Galaxy. We assume that Rdiff = 10 ^ 20 
which is appropriate for high-latitude clouds fBlit z et al.l 
[1984; Magnani et al. 1985) and smaU clouds in the ISM 
(Q^naoD & Bowers 1988). Then we can use the data in 
FigurelHlto estimate the fraction of diffuse emission in the 
different regions of the galaxy. From 200 pc < Rgai < 800 
pc, Ravg ~ 6 which implies, under these assumptions, 
that 40^50% of the CO emission comes from GMCs 
and the balance comes from translucent molecular gas, 
comparable to the fraction seen in the Milky Way. A sig- 
nificant fraction (20%) of the CO emission comes from 
the inner 200 pc of the galaxy, around the AGN. Here, 
i?i3 = 8 which implies that between 65% and 85% of the 
CO emission from the nucleus arises from diffuse molec- 
ular material which is equal to the fraction over the en- 
tire galaxy. In contrast, between 25% and 60% of the 
^■^CO emission comes from diffuse gas. Assuming the 
^■^CO emission linearly traces the molecular mass of the 
galaxy, we assume Rdiff = 20 and estimate that 25% of 
the molecular mass resides in diffuse structures and the 
remainder is in bound clouds. 

The results of the Fourier analysis of the emission sug- 
gest a similar trend. In the ^^CO emission, the power 
on small {£ < 200 pc) scales comprises 20% of the total 
power, but for CO, this decreases to 10%. The Fourier 
analysis does not completely separate diffuse from bound 
emission since the diffuse component of the molecular 
ISM contributes power on small scales as well. 

Line ratios indicate that a substantial portion of the 
molecular emission arises from diffuse gas, where the 
standard CO-to-Il2 conversion may be suspect. The 
^'^CO-to-H2 conversion factor is more appropriate for es- 
timating the mass of the galaxy. Using X2 = 1, we find 
that Mco = (3.2 ± 0.4) x 10^ M^. U sing the ^^CO- 
to-Il2 conversion factor of iLada et "all fl994) , we find 
Miaco = (2.7±0.1) x 10^ Mq, roughly 30% of the Milky 



Way value l)Damdll993| ). The two mass estimates for the 
molecular gas agree within the errors in the flux mea- 
surements, and the mass measurement of the gas appears 
robust. Using ^^CO emission to measure molecular mass 
will likely provide a more robust measurement of molec- 
ular mass in star-forming molecular clouds rather than 
the diffuse gas which is presumably inert with respect to 
star- formation. 

5.2. Dust and Extinction in M64 

Both CO isotopomers confirm a large mass of molecu- 
lar gas in the inner kiloparsec of M64; and with the im- 
plied large column density of molecular material, there 
should be accompanying dust. Bendo et al. (2003) use 
ISO and SCUBA data to measure the long wavelength 
(A > 100/im) infrared and submillimeter emission from 
M64 and find an emission distribution consistent with 
two dust components at Tkin = 14 K and 28 K respec- 
tively. The total dust mass implied by their observa- 
tions is 2.7 X 10^ Mq which implies a gas- to-dust ratio 
of ^ 100 by mass, very close to Galactic l)Bohlin et alJ 
11978ft . If the molecular gas and the corresponding dust 
were distributed uniformly, the extinction through the 
disk should be ^ Ay — 7. Stella r extinction mea sure- 
ments imply Ay ~ 1.5 (NUGA, iWitt et alV l9M] im- 
plying that the dust (and gas) must have a filling frac- 
tion smaller than unity. The interferometer data show 
{TA,max) — 1-7 K for CO but it is unlikely that the gas 
is significantly colder than T^in = 10 K or warmer than 
Tdust = 14 K. If the difference is because of beam dilu- 
tion, we estimate that the typical covering fraction of the 
gas in GMCs is 15% < / < 45%. 

5.3. Star Formation in M64 

Large surface densities of molecular gas are linked to 
high star formation rates l)Kennicutill998l : IWong fc Blitd 
12002ft . In this section, we synthesize the results of the 
observations in several wavebands to determine if the 
star forming properties of the molecular gas in M64 are 
peculiar. 

Using the FIR fluxes from the IRAS Bright 
Galaxy Sainple |Soifer et al.lll989() . and the scahngs of 
IHelou et alJ l|1988|) . we derive an FIR flux of 1.9 x lO'^ 
erg cm"2 Hz'^ from ^eo = 30.2 Jy a nd gmn = 78.7 Jy . 



Scaling this with the relationship of iKennicutq l)199g 
yields a total SFR of 0.17 Mq yr~^ for Rgai < 1 kpc or 
0.05 Mq yr^^ kpc^^, nearly an o rder of magnitude lar ger 
than in the Solar neig hborhood ("Mil ler fc ScaIollT979ll . 

We synthesized a radio continuum spectral en- 
ergy dist ribution from severa l studies in t h e lit- 
erature llCondon et all 1199^ iTurner fc Hoi 119941 



llsrael fc van der HulstI Il983t iFanti et all Il973|l and 

derive 



S{iy) = (0.130 ± 0.005) Jy 



-0.61±0.03 



.4 GHz/ 

The radio continuum emission is confined to the mo lecu- 

111111991 and 



lar gas disk IjCondon et al. 1998: iTiirner 
is well within the primary beam in all the interferometric 
studies. The resulting SFR is 0.15 ±0.04 Mq yr~i using 
the for mula of . C ondon (1992). Comparing this to the 
study of lMurgia et a l. ( 2002) indicates that M64 is quite 
similar to other galaxies with similar surface densities of 
molecular gas. 
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Ha is one of the most commonly used star formation 
tracers in the galaxy, though using it requires correcting 
for extinction, which can be large. B94 use a wide field 
image and estimate the Ha flux as 4.8 x 10~^^ exp((rHa)) 
erg s~^ cm~^, where (rua) is the typical optical depth for 
Ha emission, taken by Brau n ct al. ( 1994) to be 1.3. The 
corresponding star format ion rate is 0.28 Mq yr~^ with 
the scalings of iKennicuttI (^998). The Paa image pre- 
sented in H2.2l has a total flux of 9.3 x 10~^^ exp((TpaQ)) 
erg s~^ cm ~^. For a com parable degree of extinction as 
assumed in iBraun jt^a l. (1994), (rpaa) — 0.23 with the 
extinction law of Cardelli et al. ( 1989). This corresponds 
to a tota l star formation r ate of 0.15 Mq yr~^ with the 
scaling of lKennicuttI l)1998f) . This is not unexpected since 
the image covers a smaller area than the Ha image. The 
recombination line emission appears to accurately trace 
the star formation in the galaxy. In order for the linc-of- 
sight extinction to the H II regions to be only Ay = 1.7 
while having the gas imply a mean extinction of Ay — 7 
again requires significant clumping of the molecular gas 
or else essentially all the Ha emission would be extin- 
guished by dust absorption. 

The synthesis of observations suggest a star formation 
rate of ~0.2 Mq yr~^, or 0.06 Mq yr~^ kpc~^, over the 
molecular, star-forming disk. This agrees well with the 
hig h surface den sity regi ons of other molecul e rich g&\ayi- 
ies (jWong^fc B litz 2002t lMurgia et al.ll200l . Given the 
large density of molecular material present, the star for- 
mation activity is relatively quiescent, with a molecular 
gas depletion time of r > 1.4 Gyr. Thus, the star forma- 
tion efficiency of t he molecular clouds is sim ilar to that 
of the Milky Way ijMoonev &: SolomonllT98^ . 

5.4. The GMCs as Star Forming Structures 

In t l4.2l we argued that there are self-gravitating clouds 
of molecular gas in M64, though their properties are sub- 
stantially different from the GMCs found in the disk of 
the Milky Way. These massive clouds also appear to be 
the star forming structures in M64. In Figure El we 
indicate the positions of molecular clouds from our cata- 
log on the continuum subtracted HST image discussed in 
>i2.'2\ Across the molecular disk, there is a strong corre- 
lation between GMCs and recombination emission from 
star formation. To quantify this, we convolved the Ha 
image to the resolution of the ^'^CO data and compared 
the brightness on a pixel by pixel basis. The correlation 
coefficient is 0.74 ± 0.02 and rank correlation coefhcient 
is 0.50 ± 0.03. Performing a similar analysis with the 
Paa image yields a correlation coefficient of 0.70 ± 0.02 
and 0.65 ±0.02. Star formation is clearly correlated with 
the gas in the Giant Molecular Clouds, but measuring 
the star formation rate within individual GMCs will re- 
quire high resolution FIR observations. Unfortunately, 
the angular resolution of Spitzer observations will not be 
sufficient to correlate the FIR flux distribution with the 
location of GMCs. 

Because the GMCs in M64 differ substantially from 
those of the Milky Way, they may have significantly dif- 
ferent star formation efficiencies. However, the global 
star formation efficiency is not markedly enhanced. Stud- 
ies throughout the Milky W ay (Mooncv & Solomon 1988:- 
iLis et al. 1991: SncU ct al. 2002) show that that the star 
formation efficiency of clouds is constant with respect to 
cloud mass. On a large scale, it does not appear that 



the M64 clouds are peculiar in their star formation effi- 
ciency. This is not startling since even in ULIRGs there 
is evide nce that t he sta r formation efficiency is constant 
ijGao fc Solomonl l2004() . It is puzzling why such a large 
variation in the internal properties of GMCs results in a 
strikingly constant star formation efficiency. 

5.5. Comparison to other GMCs in High Density 
Environments 

There are several published studies searching for 
molecular clouds in high density environments. In par- 
ticular, we compare the results of M64 to a CO study o f 
the Antennae Galaxy (Arp 244) bv lWilson et all ll2003f) 
and C O observations of the Galactic center bv Toka^t^l 
()2001|) . Although the studies use different decomposi- 
tion methods, it interesting to compare the results be- 
tween the studies. We plot the mass-linewidth relation- 
ship for the studies under consideration in Figure 1151 
We choose to compare the mass-linewidth relationship 
since it avoids accurately measuring deconvolved radii of 
clouds across several studies. We include the study of S87 
for comparison with these three studies in high density 
environments. The masses derived from CO luminosity 
have been rescaled to X2 = 1. 

Although all the studies use different methods for cloud 
decomposition, we dir ectly compare their d erived clouds 
to the M64 sample. iWilson et all l)2003t) decomposed 
CO emission from the Antennae galaxies using the orig- 
inal CLUMPFIND algorithm, isolating massive molec- 
ular clouds. Performing a linear fit between log(Ay) 
and log(A/co) for their data yields an index of 2.0 ±0.1, 
consistent with the M64 results and markedly different 
from the Milky Way values. The offset in luminosity is 
nearly an order o f magnitude larger than the M64 data. 
'Oka e t al] l)2001|) decomposed the CO emission from the 
galactic center using the methods of S87. The index of 
their relationship is (2.7 ± 0.2) and the offset is roughly 
an order of magnitude lower than the M64 clouds. In all 
three cases, the clouds follow a different mass-linewidth 
relationship than do Milky Way disk clouds. 

For 8 resolved clouds in the Antenna study, the virial 
masses agreed with the luminous mass estimates (for 
X2 = 1.5) and the average surface density of the clouds is 
100 Mq pc~^, comparable to the clouds of S87. On the 
opposite extreme, the mean surface density of the molec- 
ular clouds in the Galactic center is 500 Mq pc~^, higher 
than M64 clouds by a factor of 2. These differences are 
particularly interesting since they suggest that the high 
surface densities observed in M64 are not intrinsic to high 
mass GMCs but likely depend on environment. One envi- 
ronmental factor which may determine the properties of 
GMCs is the external hydrostatic pressure. GMCs can- 
not exist as discrete entities for much longer than a sound 
crossing-time unless they are self-gravitating and gravity 
is balanced by the internal pressure of the cloud. In en- 
vironments that have high external pressures relative to 
the Milky Way, the internal pressures of clouds should 
be proportionally larger as well or else the pressure fluc- 
tuations in the diffuse ISM will destroy the clouds. To 
test the influence of environment on CMC properties, we 
plot the internal pressure of GMCs as a function of the 
external hydrostatic pressure for several galaxies (Figure 
For M64 and the disk of the Milky Way, we esti- 
mate the pressure using methods described in ^ For 
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Fig. 14. — - Locations of ^^CO molecular clouds overlaid on Ha grayscale image. Resolved clouds are indicated with virhite stars and 
circles with of radius Re and unresolved clouds indicated with black staxs. The border of the WFPC2 image is shown as a gray dashed 
line. Across the molecular disk, there is a correlation between GMCs and recombination emission, presumably from star formation. 




1 n" I 

1 10 100 

iV [km/s] 

Fig. 15. — The mass-linewidth relationship for four popula- 
tions of GMCs. The GMCs from M64 are compared with those 
found in the Galacti c center I Oka ct al. 200lJ) the Antennae G alaxy 
IWilson et al.l2003l) and in the Galactic disk ISolomon et al.>1987l 
S87). All studies in high density systems have a similar mass- 
linewidth relationship M oc AV^, which is dramatically different 
from the Milky Way relationship. The three high density cloud 
populations have significantly different offsets in the relationship. 
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Fig. 16. — A comparison of internal and external pressures for 
Giant Molecular Clouds in a variety of systems. Points represent 
pressure estimates in clouds for M64 (p resent work). Milky Way 
(S87), the Antennae galaxy IWilson et al.i2003iV and the Galactic 
center (Oka ct al. 200^. Internal pressures are estimated from 
cloud properties. External pressures are quoted values from X-ray 
studies (Antennae, Galactic center) or estimated from the mass 
contributions of galactic components (Milky Way, M64). For most 
clouds, the internal pressure scales with the external pressure and 
is consistently above the derived value. 



M64, we assume that 1/4 of the molecular mass is dif- 
fuse and therefore contributes to the external pressure 
with X2 — 0.25. We have included a point represent- 
ing the mean surface density for GMCs estimated from 



the Galactic center data ijOka et alJl200H) at an external 
pressure of Pext/fc = 5 X 10^ K cm-3 (ISpers-el fc Blitz! 
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119921) . We pl ot a point for the mean surface density of 
the sample of iWilson et alJ l|200,'^ using their estimate 
of the pressure, Pint/k = 6 x 10^ K cm~^. Finally, we 
include the cloud s of S87 and derive pressur e values us- 
ing the forrniolajjf Blitz & Rosolowskv f2 004j) . the stellar 
de nsity o f IDehnen fc Binnev. (,1998^) and the gas profile 
of IDamel ll199.'^ . We have averaged the internal pres- 
sures of the clouds in bins of APext/k = 2 x 10'^ K cm~'^ 
to highlight the significant scaling of internal pressure 
with external pressure within the S87 sample. This scal- 
ing is present using both luminous and dynamical mass 
estimates of the GMC in the Milky Way and implies a 
systematic variation of GMC column density with Galac- 
tic environment. There appears to be some offset be- 
tween these populations, but this may be due to dif- 
ferent methods for evaluating the pressure. The point 
representing the mean value for the Antennae galaxies 
is notably discrepant from the trend. This may be be- 
cause the objects cataloged by Wilson et al. (200j) are 
physically distinct from those found in the other sam- 
ples, perhaps because of the large spatial resolution of 
their study {Obeam 380 pc). Omitting the point for 
the Antennae, we find a significant scaling with pressure: 
Pint oc P°^l^^°-°^. This scahng suggests that GMCs in 
molecule-rich, high-pressure environments are character- 
istically denser than those found in the Solar neighbor- 
hood. Since GMCs in M64 also show scaling of surface 
density with molecular mass, it is unclear as to whether 
the external pressure sets the surface densities of the 
cloud or whether high pressure environments form high 
mass clouds and the scaling with surface density is set 
by the cloud mass. 

6. SUMMARY AND CONCLUSIONS 

We have presented BIMA and FCRAO observations 
of the molecule rich galaxy M64 (NGC 4826) in the 
{J — 1 ^ 0) transition of three tracers: CO, ^^CO and 
HCN. The spatial resolution of the observations projects 
to 75 pc at the distance of M64, making it feasible to 
marginally resolve large molecular clouds in the galaxy. 
We decompose this emission into clouds that are high 
mass analogues of Giant Molecular Clouds in the Local 
Group. The properties of these clouds are significantly 
different from those seen locally and are summarized in 
Table 0] We report the following conclusions: 

1. The Hue ratios M^co/W^iaco = Ri3 and 
M^hcn/W^co vary over the face of the galaxy. The large 
values of i?i3 seen in the nuclear regions of the galaxy 
{Rgai < 200 pc) and beyond Rgai = 800 pc imply sub- 
stantial emission from diffuse molecular gas. Between 
these two regions, R13 = 6, similar to the Milky Way 
disk. The ratio Whcn/W^co declines monotonically from 
0.09 in the middle of the galaxy to 0.012 at Rgai = 800 pc. 

2. A Fourier analysis of the emission distributions 
shows that the value of R13 = 4.1 on small scales {i < 200 
pc) and i?i3 = 6.5 on large scales. This implies that the 
dumpiness seen in the ^'^CO image of the galaxy repre- 
sents real cloud structures with values of R13 comparable 
to Milky Way GMCs. 

3. We decomposed the ^^CO emission into clouds and 
calculated their properties. We find several clouds with 
M > 10^ Mq larger than any cloud seen in the disk 
of the Milky Way. The clouds show L cx AF2.2±o.4 



which is a much smaller index than is seen for clouds 
in the Milky Way disk where L (x AV^. Similarly, 
we find AV cx compared to the Milky Way 

clouds where AV oc i?e'^. Despite these differences, 
virial estimates of their masses agree with those de- 
rived from the ^^CO emi ssion for th e Gala ctic ^^CO- 
to-Ho conversion factor of lLada et ail llT99l . The CO 
molecular clouds imply a COto-H2 conversion factor of 
(3.2 ± 0.2) X 10^° cm-^K km 

4. High mass GMCs in M64 are significantly smaller 
and denser than expected from extrapolating Local 
Group trends. The mean surface density of the M64 
clouds is 250 Mq pc~^ and the surface density of the 
clouds varies with cloud mass: I]h2 oc m"-^^"-^, com- 
pared to Sh2 oc const, in the Local Group. Clouds with 
masses comparable to those seen in the Milky Way also 
have comparable linewidths and radii, suggesting that 
the clouds may be similar for M ~ 10^ Mq. 

5. Diffuse gas (i.e. not bound in GMCs) contributes 
roughly 2/3 of the CO emission seen from the galaxy, 
but it only accounts for 1/4 the molecular mass in the 
galaxy. Nonetheless, empirically derived conversion fac- 
tors for CO and ^^CO emission give a similar value of 
molecular mass for the entire galaxy: 2.7 x 10^ Mq (for 
D = 4.1 Mpc) for X2 = 1 

6. Recombination line emission, FIR and radiocontin- 
uum measurements of the star formation rate M64 yield 
^SFR = 0.06 Mq yr-i kpc'^ or 0.2 Mq yr'^ for the cen- 
tral molecular disk. The molecular gas depletion time is 
^-^1.4 Gyr, comparable to similar galaxies. 

7. The GMCs found in M64 are star forming struc- 
tures, judging from their correlation with Ha and Paa 
emission in the galaxy. Despite the implied large values 
of their internal densities and pressures, the star forma- 
tion efficiency is not dramatically enhanced. The GMCs 
in M64 share many qualities with GMCs found in the 
center of the Milky Way. 

8. The GMCs in M64 share many properties with pop- 
ulations of clouds found in the inner Milky Way and the 
Antennae galaxy. GMCs in all three galaxies have a well- 
defined mass-linewidth relationship with an index ~ 2, 
markedly different from the clouds in the inner disk of 
the Milky Way. We the internal pressure of GMCs scales 
with the external pressure of the ISM: P„t oc P°^J^^° °^ 
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TABLE 4 

Summary of the GMC properties 





Milky Way 


M64 


Maximum Mass [Mq] 


6 X 10^ 


2 X 10^ 


Size— Linewidth Scaling 


Ay (X 


AV oc Rl ° 


Luminosity-Linewidth Scaling 


L oc 


Lcc AV^-^ 


Mass— Surface Density Scaling 


Sh2 oc a/0-0 


Sh2 oc MO-7 


{Eh2> [Mq pC-2] 


100 


250 


{Pi„t) [K cm-3] 


10^ 


lO'' 


SFR [ Mq yr-l kpc-^] 


0.005 


0.06 


Mol. Gas Depletion Time [Gyr] 


0.5 


1.4 



APPENDIX 

THE MODIFIED CLUMPFIND ALGORITHM 

In this paper, we use a modified CLUMPFIND l)Williams et al.llT99l WGB94) algorithm since we seek maximum 
flexibility in identifying the geometry of the molecular emission. The dataset we are analyzing in this paper is 
significantly different from that of WGB94, who used single dish data that were spatially Nyquist sampled. In 
contrast, the interferometer data has several pixels across a beam and is often at lower signal-to- noise than the data 
examined by WGB94. To optimize the CLUMPFIND method for our data, we adopt two significant changes. First, 
we determine which local maxima to include in the analysis before dividing the emission into clouds corresponding to 
those peaks. This allows a more careful analysis in the relatively low signal-to-noise domain. The other change we 
adopt is to change the path length metric so that path segments in velocity are a factor of Wy = Ofwhm /'^^x longer 
than path segments in velocity, is the number of spatial pixels per beam half-width, roughly 4 in these data. This 
change accounts for the oversampling in the spatial direction of our data. A graphic representation of the distance 
metric adopted in the datacube is shown in Figure IXT1 

Accurate selection of significant local maxima is the key to implementing our algorithm. In WGB94, significant local 
maxima are determined by contouring the data in units of AT/arms = 2, descending through the contour levels, and 
marking any new, isolated island at a given contour level as a new maximum. However, a noise fluctuation can cause 
an artificial local maximum and a false clump and the choice AT/arms = 2 is a compromise between discriminating 
against noise spikes (large contouring interval) and sensitivity to emission structure (small contour level) . To optimize 
the algorithm, we separate these two criteria: adopting a small contouring level to follow the emission structure well 
while keeping track of the levels at which isolated regions merge in order to later determine wheth er certa i n "clo uds" 
represent noise spikes. We adopt a similar method to (and heavily inspired by) the work of Bru nt et al.l l)2003tl . In 
their method, two lo cal max imum arc considered distinct if the saddle point connecting them is at least ATuniq / o'rms 
below both max ima. iBrunt e t al. (2003) use ATuniq / o'rms = 4 but we find a value of 3.0 is appropriate for our data 
(see Figure IX2|) . In our method, we contour the data in units of 0. Scrims, and record isolated islands that appear in 
each contour level as well as the level at which pairs of maxima merge. If a maximum merges with another, higher 
maximum and the saddle point between the maxima is less than ATuniq/ Ormsi it is considered to be the same cloud 
and the smaller local maximum is removed from the analysis. This selects against noise spikes since they appear as 
isolated maxima in a single resolution clement of the datacube whereas cloud structures appear in multiple elements. 
Once all connected local maxima are tested for the presence of noise, the remaining local maxima are used to partition 
the data set. 

The data set is partitioned by stepping downward from the maximum value of Ta in units of Q.bdrms- Pixels added 
at each level are assigned to the local maximum that has the shortest path through the datacube connecting the pixel 
to that maximum. Paths are constrained to run through pixels which contain emission, and the length of a path 
between two pixels is defined as the minimum length over all possible paths that lie within the bounds of the emission. 
Steps in velocity are considered to be a factor of Wy = Ofwhm /'^Ax longer than steps in position space. If pixel is 
connected to two maxima by paths of equal length, then the pixel is assigned to the closest pixel using a Euclidean 
distance metric: 

d(xi,X2) = ^/{xi - X2f + {yi ~ y2Y + wl{vx - ua)^. (Al) 

Here, Wy is the same factor that accounts for oversampling in the spatial direction. If the pixels have equivalent path 
and Euclidean distances to two separate maxima, the pixel is assigned randomly to one. 

To evaluate the optimal parameters in running the algorithm, we performed Monte Carlo simulations on controlled 
simulations. The simulated data consisted of two clouds with Gaussian profiles in position and velocity space. Inter- 
ferometer noise was added to some trials to mimic the observation conditions. Ideally, the algorithm will consistently 
separate the blended Gaussians into (only) two clouds. By varying the signal-to-noise, Gaussian separation and al- 
gorith m pa rameters, we explored the limitations of the algorithm and the ideal parameters for decompositions. In 
Figure fX2l we display the results of two such experiments. First, we fixed the cloud separation at 2>(Jx and varied 
the peak signal-to-noise to determine the limitations imposed by noise on the decomposition. For ATuniq/ Orms = 3, 
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Fig. A1. — Distances from the central shaded pixels to adjacent pixels in the data cube. The number Wy accounts for the oversampling 
in the spatial direction by making different velocity planes correspondingly "farther" from the central pixel. Within a given plane, only the 
four pixels immediately adjacent to the central pixel are considered neighbors and, consequently, one distance unit away. This definition 
results in a "taxicab metric" in a given plane. 




5 10 15 20 25 1 2 3 4 

Fig. A2. — Benchmarking results for the modified CLUMPFIND algorithm. In the left-hand panel we plot the number of clouds the 
algorithm that result from decomposing a trial data set as a function of signal to noise. The modified algorithm produces good results even 
down to a peak signal-to-noise value of 5. In contrast, the CLUMPFIND algorithm displays large sensitivity to noise in the low signal- 
to-noise regime. In the right-hand panel, we plot the number of clouds that result from decomposition as a function of cloud separation 
for Tp^ak/'^rms = 8. The trials are conducted using a range of values for the algorithm parameter l:^T^niq/'^rms, which determines the 
algorithm's sensitivity to noise. The trials show that /yT^niq/ arms = 3 is reproduces the model well while not requiring extremely high 
signal-to- noise in the data and that for this choice, the algorithm separates the two model clouds for Ax'/cr^ > 3. 



we find that our modified CLUMPFIND algorithm consistently recovers two clouds while the original CLUMPFIND 
finds a significantly larger number owing to the presence of no ise. The modified algorithm breaks down when the 
signal-to-noise ratio is less than 5. In the second panel of Figure IX2I we compare the number of clouds resulting from 
the decomposition as a function of cloud separation and ISTuniq/ (^rms- This experiment explores the ability of the 
algorithm to separate blended clouds and motivates our selection of ^Tuniql <^rms — 3. We choose a signal-to- noise 
ratio of 8 for the trials. Selecting ^Tuniql <^rms < 3 results in noise spikes causing a spurious decomposition of the 
clouds. Selecting /^Tuniq/ Crms > 3 does not significantly change the results of the decomposition. The algorithm 
successfully separates clouds with a separation between them of Ax /ax ^ 3. 

We have also checked the influence on choice of contour levels and find that a contouring the data in levels smaller 
0.5 (Jrrns docs not change the outcome of the decomposition. We use 0. Scrims in our analysis of the data. In addition, 
we also check the ratio of fluxes in the resulting clouds for trial clouds with unequal intensities; and we find that the 
algorithm partitions the clouds reasonably. 
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DERIVATION OF CLOUD PROPERTIES 

Using intensity- weighted moments, we have calculated the macroscopic properties of the clouds in M64. The position 
of the cloud is taken as the centroid position, weighted by the antenna temperature in each pixel Ti, e.g., xq = 
(Si '^i^i)/ (Si -^i)- The linewidth and radius of the cloud are calculated by taking the second moment of the emission 
in velocity and position respectively: 



cr? ^J-jT [(^ ~ ^o)^ + {y- yof 



^1 = v=r777 'Yt,{v- vo 



(Bl) 
(B2) 



To transform the moments into cloud radii and linewidths, we first corrected the derived moments for the effects of 
clipping, since the lowest contour level is at Ta — 2(Trms rather than zero. As such, the derived moments underestimate 
the true cr^ and ct^ of the cloud. In order to compensate, we assumed the clouds have a Gaussian profile in position- 
position-velocity space. Since the clouds are only marginally resolved at best, the assumption of a Gaussian profile 
seems to be reasonable. We correct the moments by the ratio of the value of th e moment extrapolated to a co mmon 
reference value to the value of the moment for emission clipped at Ta = 2arm.s l|Oka et al.ll200lt iBolatto et al.ir20031 
The correction factor depends solely upon the peak-to-edge ratio, P = Tmax/TcUp where Tmax is the maximum 
antenna temperature in the emission region and TcUp is the clipping level, equal to 2arms for this study. In this case, 
the correction factor, f{P), is 



f{P) = 



J^"" exp(-a;V2) 



Jo 



■ exp(— x^/2) 



x^exp(-xV2) 
/p """^ x'^ exp(— a;2/2) 



1 -1 



where 



(B3) 



= A/21n(P) and Xref = 2 



The clouds in iSolomon et al.l l)1987|) typically have P — 8 implying, for that sample, Xmax = ^2 ln(8) = 2.04. To 
match the clipping levels used in that study, we choose the reference value for the extrapolation Xref — 2. For a typical 
cloud which has P = 3.5, f{P) = 1.4 resulting in a correction of 20% to the size and linewidth of the cloud. 

The linewidth is calculated by scaling the corrected dispersion to a FWHM using the Gaussian linewidth correction: 
AV = <7y ■ ^ fiP)8 lii(2). The cloud radius is derived by correcting (7^ for clipping and beam convolution and then 
applying the scaling from ct^ used bv ISolomon et all l)198 

3.4 



in their study of Milky Way clouds: 



Re 



/27r 



2(7,2 _ 

beam 



(B4) 



Our simple model to correct for beam convolution breaks down when ar\/f{P) < (Jbeain, and so we only consider 
a cloud resolved if Re > 40 pc. Similarly, we only consider a cloud to be resolved in the velocity direction if AV is 
larger than the channel width. 

The luminosity of the cloud is also reduced by clipping the cloud at 2arms ■ To account for this clipping, we correct 
the integrated intensity value by a factor of 



9iP) 



JO 



/p'''"°"a:2exp(-a;2/2) 



(B5) 



For clouds with P = 3.5, the correction factor g{P) — 1.4 for Xref = 2. The dominant source of uncertainty comes 
from the extrapolation of cloud properties using the correction factors f{P) and g{P), though most of our conclusions 
are independent of applying this correction. In this paper, we have assigned the derived quantities an uncertainty 
equal to half the value of this correction. 
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